Weak localization effects in granular metals 



o 
o 

o 



C. Biagini\ T. Caneva^, V. Tognetti^-^'^ & A. A. Varlamov^ 
^ INFM, Unita di "Tor Vergata", Viale del Pohtecnico 1, 00133 Roma, Italy 
^ Dip. di Fisica, Via G. Sansone 1, 50019 Sesto F.no (Firenze), Italy 
^ INFN, Sezione di Firenze, Via G. Sansone 1, 50019 Sesto F.no (Firenze), 
* INFM, Unita di Firenze, Via G. Sansone 1, 50019 Sesto F.no (Firenze), Italy 
5 GOHERENTIA-INFM, GNR, Viale del Politecnico 1, 00133, Roma, Italy 
(Dated: February 2, 2008) 

The weak localization correction to the conductivity of a granular metal is calculated using the di- 
agrammatic technique in the reciprocal grain lattice representation. The properties of this correction 
are very similar to the corresponding one in disordered metal, with the replacement of the electron 
mean free path I by the grain diameter d and the dimensionless conductance g by the tunnelling di- 
mensionless conductance (Jt- In particular, we demonstrate that at zero temperature no conducting 
phase can exist for dimensions D < 2. We also analyze the WL correction to magnetoconductivity 
in the weak field limit. 
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Recently, the properties of granular materials have at- 
tracted special attention (HI, j2|, 1^ [J] . The quantization of 
the electron spectrum in small grains requires to revise 
the basic idea of the quasiparticle spectrum continuity, 
assumed in the description of most properties of metallic 
and superconducting systems. The appearance of a new 
energy scale, the mean level spacing S, results in unusual 
superconducting properties of such systems, possibility 
to observe there specific quantum phase transitions, etc. 
In particular, it turns out that the interplay between the 
intra-grain diffusion and inter-grain tunnelling of elec- 
trons makes the metal-insulator transition in such " quan- 
tum metal" very peculiar. In this Letter, we intend to 
discuss the specifics of the weak localization corrections 
in granular systems. 

The elastic electron relaxation rate in granular metal 
consists of three contributions 
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where 



imp 



is the mean scattering time of electrons with 



impurities, Et — vp/d is the Thouless energy, d is the 
characteristic grain size, vp the intra-grain Fermi veloc- 
ity). The last term is the electron inter-grain tunnelling 
rate: T ~ gr^, where gr — Stt (t/S) is the tunnelling di- 
mensionless conductance and t is the tunnelling energy. 

The character of electron motion in a metal is con- 
veniently classified as a function of its mean free path 



£0. The diffusion length 



y/Vn/T (where V„ 



VpTimpI D is the diffusion coefhcient of metal and D is 
the space dimensionality) separates the regions of bal- 
listic :s> ir) and diffusive {£ <C It) electron motion. 
When £ h/pp the metal-insulator transition in 3D 
case takes place. In the case of a granular metal, one 
can expect that the standard WL theory, describing the 
precursor effects of this transition has to be mod- 

ified in two ways. First, the diffusion coefficient V here, 
at least in some interval of parameters, has to be de- 
termined by the tunnelling time F""'^ instead of elastic 
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FIG. 1: A schematic representation of the two possible sce- 
narios for the electron motion in granular metals as a function 
of the ratio between d, the grain diameter, £t, the thermal 
length, and £, the mean free path. 



scattering time Timp and therefore we will have a dif- 
ferent tunnelling thermal length £t — ^/DjT . Second, 
the characteristic grain size d must appear in the theory 
side by side with the diffusion length Ip. Two differ- 
ent situations are possible. When d ^ £t, the diffusive 
part of the axis is divided to the regions of the normal 
intra-grain diffusion {h/pp ^ £ ^ d) with the metal- 
lic diffusion coefficient 2?„ and of the tunnel inter-grain 
electron diffusion {d ^ t ^ £t) characterized by the dif- 
fusion coefficient V — (^(x{t)^'^ /tj ^ Fd^ (see Fig. 

^). In the opposite limit, when d ^ £t, only the intra- 
grain diffusion is possible, but, in its turn, the domain 
of ballistic regime contains two regions: the intra-grain 
one with ballistic time Timp {£t ^ ( ^ d) and the region 
of inter-grain electron motion with the ballistic time F~^ 
(see Fig. Ut). 

We will discuss here the most interesting case of low 
temperatures, d <^ Ip- We will also assume that the 
electron motion inside a single grain is ballistic {tbuik ^ d, 
£ ^ d): this means that before tunnelling to the neighbor 
grain the electron is reflected many times on the grain 



boundaries (F ^ Ep). As a consequence, Tg, 



E^ 
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and the dimensionless conductance of a single grain g = 
Et/5 ^ gr- The conductance of the entire system is 
given by g^^ + g^^ ~ 9t^j what is equivalent to say that 
the drop of applied electric potential occurs only inside 
the tunnel barrier. 

We assume that the grains are almost identical, with 
average diameter d, and form a regular lattice with lattice 
constant equal to the same d. The coordinate of the grain 
center will be labelled by the lattice variable R^. The 
Hamiltonian of the system can be written as 



H 



i,p 



h.c. 



where Ci^p (c|p) is the annihilation (creation) operator 
of an electron in grain i with intra-grain momentum p. 
The tunnelling energy tfj^ will be taken equal for all 
bonds between nearest neighbor grains and independent 
on the intra-grain momentum, i?'^ = t. Performing the 
Fourier transform with respect to such one can write 
the Hamiltonian in the representation of both intra- and 
inter-grain momenta (double momentum representation) : 

K,p 



K p^p' 



ck,pCK,p' 



h.c. 
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Here K is the quasi-momentum belonging to the recipro- 
cal to Ri lattice and it varies in the first Brillouin zone. 
The lattice structure factor 7k Z^^ J2^=i e*^ '^*'; 
are the vectors connecting the center of selected grain 
with the nearest neighbor sites, Z is the coordination 
number. For a simple cubic lattice, the vectors have 
one component equal to ±d and all the others equal to 
zero. For sake of simplicity, we will restrict to this case; 
the extension to generic lattices is straightforward. In 
simple cubic lattices, 7k = j)^a=i'^'~'^(^a'^)- From 
Eq. we can define the single electron Green's func- 
tion in the double momentum representation as: 



GK(P,£n) = 
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iSn - Cp - Zt{l -7k)' 

ttT (2n -I- 1) as 



with e„ = e„ + (2Tei) signe„ and e„ 
fermionic Matsubara's frequency. 

Recalling that the electric field is negligible inside the 
grains and differs from zero only inside the barriers, in 
presence of the vector potential A one can write the a— th 
component of the electrical current operator in the imag- 
inary time T as 
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FIG. 2: Dyson equation for the diffusion vertex part in a 
granular metal in the direct space. The thick lines are the 
double Fourier transform of the fermionic Green's functions 
Eq. @, the dashed lines represent single impurity scattering 
and the initial and final grains are reported as roman indices 
i, j etc. The Cooperon is obtained from the diffuson via time- 
reversal transformation on a fermionic Green's function _6j. 



The linear response function, expressed as the second 
derivative of the partition function, is given by 
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Z[0] M„(r)Ma'(0) 



Ila,o.'iT)-c'd'5^^,S{T){HT)^, (5) 



where the current-current correlation function is ex- 
pressed via the current operator as 



n^^a' (r) = -(r,J„(r) (0) 



The thermal average (...)o shall be performed with the 
diagonal Hamiltonian, the first line of Eq. (0) . We are in- 
terested in the diagonal components of conductivity ten- 
sor: 



Iia,a{uJ^) = 2e2d2 lil^^sin^ (ifad) 
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^ r^^GK(p,£„+.)GK(p',e„). (6) 

e„ p,p' 

where En+v — £n + and uj^ — 2t:Tv is the bosonic 
Matsubara's frequency. We can formulate the following 
rules of diagrammatic technique in the double momen- 
tum representation: 1. at each external vertex attach 
a factor = etdsm{Kad)\ 2. at each straight line at- 
tach a single electron Green's function Gk (p, £n)] 3. sum 
over all internal momenta and Matsubara's frequencies; 
4. impose energy and lattice momentum conservation at 
each vertex. Now the impurity averaging of Eq. lO can 
be performed. The Cooperon vertex corresponding 
to the granular metal can be obtained from the Dyson's 
equation reported in Fig. Other possibility is to 

renormalize the standard intra-grain Cooperon by means 
of introduction in the corresponding diagrams of a self- 
energy correction appearing due to tunnelling, as it is 
done in Ref. 0. Both approaches in the assumptions 
made above turn out to be completely equivalent and 
lead to the expression: 
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FIG. 3: WL correction to the conductivity in the double mo- 
mentum representation. The solid lines are single electron 
temperature Green's functions Gk(p, £n); the external ve- 
locity vertices are Va = etdsm{Kc,d); the shaded box is the 
Cooperon Cq [lou)- 



In the latter expression appears the exact value for the 
tunnelling rate: V = ZqtS. In the expression V = 



2r(l-7Q)/|Q| 



TcP one can also recognize 



the effective " tunnelling diffusion constant" . 

Now one can directly calculate the diagram reported 
in Fig. In our assumptions, see Eq. the Green's 
function can be considered independent on K when inte- 
grating over p, because its behavior is completely deter- 
mined by the pole due to the impurity (or grain bound- 
aries) scattering, related to rj^, just as in the case of 
the diffusive limit for a bulk system. Performing the fre- 
quency summation, the p and p' integration and the sum 
over the lattice momentum K one finds 
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At this point, we can find the WL correction to the con- 
ductivity as 
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In the case of a bulk granular system [D ~ 3) the 
sum converges and the correction is finite. The metal- 
insulator phase transition can be observed at a critical 
value Zg"^ - 0(1). 

In the case of granular film or wire [D < 2) the WL 
correction Eq. ((Jj) diverges at small Q — >0. This fact 
indicates that at zero temperature these systems cannot 
exhibit metallic properties for any value of the dimen- 
sionless conductance gx- At finite temperatures, a natu- 
ral cut-off of Eq. |7J is provided by the phase-breaking 



6]. It is related to the 



rate 7^ = {TT^y^ = (d/L^) 
phase-coherence length L,^ — yyVr^, which tends to in- 
finity when T — !■ 0. In result 
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Eqs. lO permit to define the localization length at 



which the correction becomes of the order of cj^j) 
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It is worth to stress that Eq. {T)) does not contain t^j^ 
but only lattice parameters as gx and the coordination 
number Z: the intra-grain dynamics, as expected, sim- 
ply drops out of the calculation. The conductivity and 
its corrections are related to the diffusion on the grain 
lattice, and the mechanism of the momentum randomiza- 
tion between different grains is not crucial: the electron 
dynamics at low temperature can be thought as that of 
a random walker on a lattice. This picture is fully con- 
sistent with the existing WL theory ^ , and it is in 
agreement with previous experimental findings in granu- 
lar metals 

More intriguing are the properties of the WL correc- 
tion to the magnetoresistance of granular metal. It is 
well known that quantum corrections to conductivity in 
disordered metal are very sensitive to the magnetic field: 
in fact, its presence disturbs the phase coherence of elec- 
trons moving along the self-intersecting trajectories, sup- 
pressing the WL correction and leading to the appearance 
of the anomalous negative magnetoresistance I3. In the 
following, we will show how such a correction manifests 
itself in the case of granular metal. 

To calculate the WL contribution to magnetoresis- 
tance, it is necessary to rewrite Eq. ^ in the direct 
space: 
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Zgr 



which is independent on i, depending only on the inter- 
grain spacing d. a represents the bond along the direction 



of the current. 



Here Cq^ 



-lUJ 



2r (1 — 7q). Let us notice that this form underlines the 
fact that transport is due only to the potential drop inside 
the barrier separating two grains i and i + a. In the 
presence of a magnetic field, the Cooperon wavefunction 
is given by the solution of the equation 



(4^)(l-7Q+2eA)^^ {r)^Ei;, (r) : 



(9) 



Moreover, also the intra-grain Cooperon will be renor- 
malized by the presence of the magnetic field, acquiring 

a mass term equal toQ Sq (H) = | ^-^j Et, where 

d> = Hd^ is the magnetic flux threaded through a single 
spherical grain and (j)Q = n/e is the flux quantum. When 

the field satisfies the inequahty d <^ £h — {^H) ' , or 
(f) <^ 4>Q, we have 

C^,{V,V',L0) (10) 
V'».Q||Q^»(r)^|,Q||Q^n(l-0 
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where fic = 4D/£|j is the Cooperon cyclotron energy. Q\\ 
is the momentum along the magnetic field and (n, Q _\_) 
are the quantum numbers of the Landau basis. 

The most interesting case is the two-dimensional ge- 
ometry, with the magnetic field applied across the plane 
of the sample^ for which the WL correction is: 
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Since the sum over Landau levels is evidently divergent 
at the upper limit we introduced the cut-off parame- 
ter riniax = TT^r/Slc S> 1, at which the cyclotron fre- 
quency becomes of the order of the zone edge, rimaxf^c ^ 
V {-K/d)^ . The final result reads 
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where we expressed the magnetic field in terms of the 
magnetic flux. 
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where ■0(2;) is the digamma function. The WL correction 
to the magnetoresistance 6(j (H) = Sa^^ {H)-6a^'^ (0) 
is obtained as 



5a (H) 
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One more energy scale shows up in Eq. (|12|l with respect 
to the bulk case, namely the Thoulcss energy. In the limit 
of very weak fields. 



bo <C yJl/^ErT^, this energy scale 
is not observable in the magnetoresistance: the leading 
singular correction reduces to 



8<7 (H) 
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which corresponds to the anomalous magnetoresistance 
of the standard theory. The granular behavior deviates 
from the bulk one at fields such that 0/0o ^ ^/Et 
where the intra-grain term starts to dominate in the sec- 
ond digamma function in J^{H,^^): for ^/l/4ETT^ <C 
0/00 <^ "ftp, the magnetoresistance correction acquires 
the logarithmic form 

(To nZgr ^ V F 0o y 

Larger fields are out of the range of our approximated 
approach, in which the intra-grain Landau levels, with 



the spectrum ojc{n+l/2) with cUc = 4'D„/£|^, start to 
significantly contribute to the Cooperon wavefunction. 
Let us notice finally that the ultraviolet cut-off rij^ax ^ 
Et/lOc — r/ric remains the same. 

In summary, we developed a diagrammatic technique 
in a double-momentum representation for transport in 
granular metals. Using this technique, the weak local- 
ization corrections to the conductivity arise in a natu- 
ral way and an explicit calculation shows the same low- 
temperature behavior as in bulk metals, but with the 
diffusion constant I?„ replaced by the effective tunnelling 
diffusion constant T) = Yd? and the mean free path £ by 
the average grain diameter d. Our result agrees with Eq. 
(13) of Ref. in the Q ^ limit; however, our tech- 
nique underlines the presence of the grain lattice, repre- 
sented by the cosine factor in Eq. ^ , reminiscent of the 
lattice structure factor 7q. We also give an estimate of 
the magnetoresistance correction for very weak fields. 
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